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Abstract 

The asymptotic behavior of geometry near the boundary of maximal 
Cauchy development is studied using a perturbative method, which at 
the zeroth order reduces Einstein's equations to an exactly solvable set 
of equations — Einstein's equations with all "space" derivatives dropped. 
The perturbative equations are solved to an arbitrarily-high order for the 
cosmological spacetimes admitting constant-mean-curvature foliation that 
ends in a crushing singularity, i..e., whose mean curvature blows up some- 
where. Using a "new" set of dynamical variables (generalized Kasner 
variables) restrictions on the initial data are found that make the zeroth- 
order term in the expansion asymptotically dominant when approach- 
ing the crushing singularity. The results obtained are in agreement with 
the first order results of Belinskii, Lifshitz and Khalatnikov on general 
velocity-dominated cosmological singularities, and, in addition, provide 
clearer geometrical formulation. 

1 Introduction 

The set of all solutions of Einstein's equations is too large from both the math- 
ematical and physical standpoint. Many attempts have been made to formulate 
the conditions that would restrict the solutions of Einstein's equations to a class 
that is: 

• large enough to include all "reasonable" physical situations, 

• small enough to exclude those deemed physically "pathological" , 

• mathematically as simple as possible. 
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The most promising candidate to date is (a generic subset of) the class of 
maximally-extended globally-hyperbolic spacetimes. The most serious problem — 
the problem that initiated the Strong Cosmic Censorship conjecture — for that 
class stems from the Hawking-Penrose-Geroch singularity theorems, which indi- 
cate that the globally-hyperbolic spacetimes, due to the possible incompleteness 
of non- curvature- singular timelike geodesies, might have to be extended beyond 
the boundary of their maximal Cauchy development. Possible incompleteness of 
curvature-singular geodesies is physically acceptable for two reasons. The first is 
that being crushed or torn apart by infinite tidal forces is believed to be a good 
reason for any observer to disappear, and the second is that there is some hope 
of rescue by a quantum modification of classical general relativity; as it is be- 
lieved that quantum effects become dominant in the regions of strong curvature, 
and might eliminate the unpleasant classical infinite-curvature singularities. 

It is, therefore, crucial to establish how often and under what conditions 
the non-curvature-singular (bad), as opposed to the curvature-singular (good), 
incomplete timelike geodesies occur near the boundary of the maximal Cauchy 
development; making the study of the asymptotic behavior of geometry near the 
Cauchy boundary crucial for the acceptability of the class of maximal globally- 
hyperbolic spacetimes in general relativity. Unfortunately, due to the complexity 
of Einstein's equations, an exact or useful description of the asymptotic behav- 
ior seems to be unattainable, at least so in the forseeable future. Some results of 
the work done by Belinski, Khalatnikov and Lifshitz hereafter BKL, indicate 
that, generically, the approach to the boundary of the maximal Cauchy devel- 
opment (the singularity in their terminology) is extremely complex and can be 
thought of as a sequence of pointwise Mixmaster-like transitions between Kasner 
epochs. A satisfactory geometrical formulation of their method is still lacking, 
and it is still controversial whether their method can give any information about 
the global structure of the singularity M. Nevertheless, an important discovery 
of BKL was that the asymptotic dynamics simplifies significantly for spacetimes 
for which, in a suitable foliation, the spatial derivative terms in the equations of 
motion can be neglected asymptotically. Since Einstein's equations thus trun- 
cated are exactly solvable, we can extract all the asymptotic properties from 
the solution of the truncated equations, the Generalized Kasner Solution; a so- 
lution which evolves pointwise like the Kasner solution. Most importantly, the 
Generalized Kasner solution is generically curvature-singular, as there is only 
one, (1,0,0), non-curvature-singular Kasner solution. 

The question I shall try to answer here, without going into technical details, 
is: Under what conditions is that approximation, usually called the Velocity- 
Dominated Approximation (VDA) H , valid in the class of globally- hyperbolic 
vacuum cosmological spacetimes, when approaching the boundary of the maxi- 
mal Cauchy development? 
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2 Generalized Kasner Variables and Velocity- 
Dominated Approximation 

Let K? and be the extrinsic curvature and spatial metric tensors induced by 
an arbitrary foliation S t on the spacetime, and let r = Tr{K) and K/ = K/ jr. 
There exists a unique triad ES\ of spacelike vector fields which puts K? and 
Qij into the following form: 

EgE T = I t- 2 ^ I , (1) 
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provided Pi ^ P2 ^ P3 an d t > 0. In that case there is a one to one map- 
ping between the old dynamical variables (if, g) and the new ones (E, T,p) with 
£ p g = 1. If we restrict ourselves to a subclass of cosmological spacetimes that 
admit constant-mean-curvature (CMC) foliation, which is unique if it exists 0, 
we obtain what I call generalized Kasner variables (GKV): a triad of spacelike 
vector fields E^ and three Kasner exponents p q (two of them independent), 
invariantly defined throughout the spacetime. Since the globally-hyperbolic 
spacetimes with surfaces of infinite mean-curvature (crushing singularities) have 
Cauchy boundaries exactly at those surfaces, we can analyse the Cauchy bound- 
aries by taking limit r — > ±oo. In spacetimes with no crushing singularities, 
however, there is no generic way to lock onto the Cauchy boundary, whose elu- 
siveness prevents us from analysing the geometry near it. So we should restrict 
our attention to the crushing singularities. 

In order to study the behavior of metric near the spacetime's Cauchy bound- 
ary, I use a perturbative VDA method previously developed by Moncrief and 
myself. It was applied successfully to T 3 x R Gowdy spacetimes, where it was 
found, by an inductive argument, that the VDA is valid to an arbitrarily high 
perturbative order ||. The method is applicable to spacetimes with no isome- 
tries as well, and, in a nutshell, consists of modifying Einstein's equations by 
multiplying the shift vector in the ADM action by an e, 

I ADM -> I AD Me = / [N H + eN*Hi] , (3) 

and then expanding the solution of the modified Einstein's equations with re- 
spect to that e, about e = 0. The introduction of e is just a device to obtain 
a well defined series of recursively solvable perturbative equations which start 
with the generalized Kasner solution as the zeroth-order term. The case we are 
interested in is e = 1. 
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After rewriting Einstein's equations in terms of the GKV, and choosing the 
time t to be some monotonous function of the mean curvature r, we find that 
in the zeroth order the GKV are time independent, that the lapse function N 
depends only on time t, and that the Hamiltonian constraint gives an additional 
restriction on (generalized) Kasner exponents, which together with the initial 
trace constraint gives just the standard Kasner relations for the zeroth order 
functions (only of "spacelike coordinates"): 

S V \ = E Pq = 1 (4) 

The zeroth order solution is just the generalized Kasner solution. 

All higher order terms, in turn, are uniquely determined by the initial data 
for the zeroth order solution and, using an inductive argument, it can be showed 
that, as t — > oo, they become negligible when compared with the zeroth order 
provided 

E x (V x E x ) = 0, (5) 

where E\ is the eigenvector with the only possible negative Kasner exponent. 
Roughly, the rate of decay is faster the higher the order. 

3 Conclusions 

The condition (^) has already been established by BKL Q as a sufficient con- 
dition for the decay of the first-order correction. Using the generalized Kasner 
variables I showed that (|J) is also sufficient condition for the decay of all higher 
orders in the expansion, asymptotically as r — * oo, and, therefore, for the VDA 
to be valid perturbatively. Even though I did not use a synchronous foliation, 
the fact that the lapse function N asymptotically approaches a spatially ho- 
mogeneous function indicates that the CMC foliation asymptotically becomes a 
synchronous foliation when approaching r = oo. 
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